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Abstract 

We consider the minimal low-energy action for gravity up to six 
derivatives which is renormalizable at the two-loop level modulo higher 
derivative corrections. Then we study the classical solutions corre- 
sponding to the Schwarzschild and the Robertson- Walker metrics. In 
the first case we find a singularity close to the gravitational radius and 
in the second case we find inflationary de Sitter solutions in absence 
of any matter or cosmological constant. 
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Introduction 

Today most of the physicist agree that the issue of quantum gravity is 
probably the most important open problem in theoretical physics (for a recent 
review see [1] and references therein). At the present moment it is not clear 
if the formulation of the quantum theory of gravitation will require going 
beyond quantum field theory like it happens in string theory or even a deeper 
modification of quantum mechanics. 

Of course this question will not be solved here, but until it is clarified one 
can adopt a much more modest attitude trying to figure up which could be 
the low-energy impact of the real theory of gravitation. It is clear that at low 
energies this theory will lead to the standard Einstein lagrangian. However, 
this lagrangian will in general be modified by higher-order terms in the sense 
of the number of derivatives of the space-time metric over the Planck mass 
scale Mp. Thus, the Einstein lagrangian with two derivatives and one factor 
Mp would be the first term, then we would have terms with four derivatives 
and one M P = 1 factor, next we would have six derivative terms with one 
Mp 2 factor and so on. The form of all these terms is dictated just by the 
number of the derivatives and general covariance. Different underlying grav- 
itational theories will predict different values for the adimensional constants 
affecting these terms, being superstring theories a typical example [2]. 

Independently of the concrete value of the constants in this low-energy 
effective lagrangian for gravitation, it is clear that one can use this derivative 
expansion up to some number of d 2 /Mp powers to predict the qualitative 
behavior of the modifications to the standard Einstein gravitation theory. 
However, in order to be completely consistent one should also take into ac- 
count the effect of the quantum corrections to this effective lagrangian. For 
example the loops obtained with the Einstein lagrangian produce corrections 
which are of the order of (d 2 /Mp) 2 , two-loop diagrams give rise to corrections 
which are of the order of (d 2 /Mp) 3 etc. 

This framework for the effective description of the low-energy behavior 
of gravity has some similitude with that proposed by Weinberg some time 
ago [3] for the description of the low-energy pion scattering. In that case the 
effective lagrangian is also the most general lagrangian compatible with the 
(approximate) chiral symmetry of the strong interactions SU(2) L x SU(2) R 
which is supposed to be spontaneously broken to the group SU (2) L+R . Pions 
are understood as the Goldstone bosons associated to this symmetry breaking 
and this is why they are (nearly) massless. The corresponding fields 7r a are 



2 



just coordinates parametrizing the coset space SU(2)i x SU(2)r/ SU(2)l + r. 
This space is a three sphere i.e. a maximally symmetric space having the 
chiral group SU(2) L x SU(2) R as isometry group. At low energies the pion 
dynamics is controlled by the effective lagrangian 

C eff = ± gab (n)d^ a d»n b + a(g ab (n)d^ a d»n b ) 2 + 

(1) 

where g a b is the metric on the sphere. Here again the quantum corrections 
to the two derivative term are of the order of 0(<9 4 ) etc. In principle the 
constants a, (3... can be obtained from an underlying theory of the strong 
interactions like QCD. Nevertheless, even in the case that this theory were 
not known or the computation of the constants were hard or impossible, the 
above lagrangian can be very usefully used by obtaining the constants from 
the experiment and using them to make new predictions. This program was 
carried out by Gasser and Leutwyler in a remarkable set of papers [4] to the 
one loop level and ever since this kind of approach to the low-energy hadron 
interactions has been followed by many particle physicists with great phe- 
nomelogical success. More recently, a similar approach has been also applied 
to the phenomenological description of the unknown symmetry breaking sec- 
tor of the standard model and in particular to the parametrization of the 
elastic scattering of the longitudinal components of the weak bosons at the 
TeV scale [5]. 

Following this very general line of thinking one could try to make an 
effective low-energy description of gravity. However, some differences appear 
when we compare pions or longitudinal components of the weak bosons with 
gravitons. First we do not have any coset space to make a Goldstone boson 
like interpretation of gravitons (see [6] for some attempt in this direction). 
Second, gravity is a gauge theory and we need to go through the standard 
procedure of quantization of these theories including gauge fixing and ghost 
terms in the lagrangian to relate properly the ten independent components 
of the metric tensor with the two physical degrees of freedom of the graviton. 
In spite of these facts, some analogies still remain with the case of the non- 
linear sigma models. For example, as it was mentioned above, the one-loop 
radiative corrections to the Einstein lagrangian have four derivatives. The 
divergent terms of these corrections were computed some time ago by 't 
Hooft and Veltman [7]. Using the background field method and dimensional 
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regularization it is possible to find: 



The important point here is that these divergences do not contribute to the 
S matrix elements (note that they vanish on shell) so that, the Einstein 
theory of gravitation considered as a standard quantum field theory is one- 
loop finite. In fact there is a simple argument to see that it must be so. 
The four-derivative local terms that in principle can be generated by the one 
loop computation are R 2 , R a/3 R a/ 3 and R a/3 ^ s R l5 af3 . Only the last one does 
not vanish on shell, but as it is well known the combination — AR^^R^ 6 a/3 + 
16R a PR a p-4R 2 is a total derivative in four dimensions related (for Euclidean 
signature) with the Euler class of the space-time manifold. Thus, modulo a 
renormalization of this topological term, the one-loop divergences of pure 
Einstein gravity vanish on shell. 

More recently, the two-loop divergent contribution to the effective action 
for gravity was computed by Goroff and Sagnotti [8]. In this huge compu- 
tation they were forced to use an special C code since standard algebraic 
manipulators were unable to work it out completely (see [9] for a previous 
disscussion about this point). The on-shell contribution to the effective la- 
grangian coming from the Goroff and Sagnotti computation is given by: 

r (2) _ 209 a/3 s p 

Leff ~ 360(4^ M 2 e K ^ K °p k «f> W 

Therefore, this remarkable computation showed that the Einstein theory of 
gravitation is not renormalizable in the common sense. However, one can 
still consider it as the first term of a derivative expansion by following the 
philosophy of the phenomenological lagrangians. The new terms that must 
be added to some given order in the number of derivatives play a double role. 
First, they have coefficients that carry the information about the underlying 
theory of gravitation and second, they are also counter terms that absorb 
the divergences. It could be possible in principle to compute the low-energy 
action of gravitation in this way up to some number of derivatives i.e. powers 
of the external momenta, in terms of some renormalized constants that in 
principle could be obtained from the experiments or observations.. Once they 
were fitted, this effective action could be used to make new predictions 



4 



Obviously, this approach has many difficult problems that must be solved. 
First, from the theoretical point of view, the practical computation is very 
hard. Note that we only know the divergent part of the one and two-loop ef- 
fective action but we even do not know the one-loop finite part contribution. 
This kind of contribution is in general non-local since gravitons are mass- 
less. Second, we have of course phenomenological problems. As the higher 
derivative terms are suppressed by powers of M P they become practically un- 
observable. The only possibility for these terms to be relevant is in physical 
contexts where the curvature is very large. Typically these situations occur 
in the regions close to the singularities appearing in the classical solutions 
of the standard theory and in fact it could be the case that the effective 
low-energy theory of gravitation completely avoids these singularities solving 
one of the old problems of the Einstein theory (we will see later that this is 
not necessarily the case). 

Motivated by this interesting possibility we have started the program 
of studying the effect on the classical solutions of the gravitation theory 
including the higher local derivative terms that appear in the low-energy 
effective action for gravity. In fact some of these terms have already been 
considered in the literature, typically the four derivative ones. However, here 
we will adopt a different attitude. As the general analysis is in principle very 
complex we want to select the minimal set of terms that must be added to 
the standard lagrangian to have a sensible low-energy effective lagrangian. 
From the Goroff and Sagnotti work we know that, at the six derivatives 
level, the only term we need to add to get an on-shell finite effective action 
is that of eq.3. Therefore, we can conclude that the addition of this term 
to the standard lagrangian produces the minimal consistent effective low- 
energy two-loop renomalizable lagrangian in the sense of the effective theories 
considered in [3] i.e. neglecting higher derivative terms. For the sake of 
simplicity this will be the only one considered in this note, so that we pass 
by the more commonly considered in the literature four derivative terms. In 
this note we will explore the effect of the six derivative term on the important 
Schwarzschild and Robertson- Walker metrics. 
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The equations of motion 

According to the above discussion we will consider the effective action: 

r e // = / d 4 x^C eff (4) 
where the effective lagrangian is given by: 

c ^ = -^ R+ w P Ra ^ w& ^ P ^ (5) 

where R a ^ 1 $ is the Riemann tensor, R is the curvature scalar (we will 
follow the notation convention of the Weinberg's book [10]), Mp is the Planck 
mass and a is an adimensional coupling constant. By applying the standard 
methods we can obtain the corresponding equation of motion: 

^(RT - \<TR) = -^(^R^sR" 5 ,^ (6) 

- ?>R a %R^ ap R^: + 6(R» 5 ap R^). AP ) 

When matter is present, we have to include a new piece in the effective 
lag as a consequence of that we must add a new term in the 

RHS of the above equation of motion which is —T^ v /2 being T^ v the matter 
energy-momentum tensor. Obviously the Einstein equations of motion are 
obtained just by taking a = 0, but for general a, the modified equations will 
produce new effects that will be studied in next section. 



The modified Schwarzschild metric: 

The most general form of the static and isotropic metric can be written 

as: 

dr 2 = B(r)dt 2 - A(r)dr 2 - r 2 (d9 2 + sin 2 6d<p 2 ) (7) 

By inserting this metric in the eq.6, we can find the equations to be satisfied 
by the A(r) and B(r) functions. For the case a — 0, these equations were 
solved by Schwarzschild in 1916 for the case of a flat spatial infinity: 

v 1 rM 2 P y ' 

A(r) = 1/B(r) 
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being M the total mass of the body producing the gravitational field. For 
general a, the task of finding the corresponding solutions is very hard. How- 
ever, following the phylosophy of the effective lagrangians, one could try to 
find at least the first coefficients of the large distance expansion of the A(r) 
and B{r) functions: 

oo r 
n=0 ' 

OO T 

n=0 ' 

where L p is the Planck lenght [L v = l/M P ) and a n and b n are adimensional 
coefficients with a — b — 1. Introducing these expansions into the equations 
of motion, we find that the first correction to the Schwarzschild solution 
appears to be of the order of (L p /r) 6 . To find these corrections, we can 
introduce the Schwarzschild solutions in the RHS of the equation of motion 
and a completely general static, isotropic metric in the LHS. With the help 
of a REDUCE code we find the following equations: 

Mp 1 A' 1 

~ 16tt VAB ~ rA 2 B ~ r^B ~ ( ' 

a 2M ^ 2352(M/MJ,) 3 _ 1080(M/M|,) 2 



My rMl' v r 9 

Ml , -B' 1 1 



16ttVA 2 B r 2 A 2 r 2 A 



+ 77t)= (11) 



a _ 2M 2A0(M/M 2 ) 3 _ 216(M /M 2 ) 2 

mJ { ~7mI )( 7$ \ 



lp I lVlp 



Ml { ^_ B> B" {B'f A'B' 

16n [ 2r 3 A 2 2r 3 AB 2r 2 AB 4r 2 AB 2 4r 2 A 2 B 1 ' 

a / -1488(M/M|,) 3 648(M/Mj,) 2 

' 77 / 



r 2 Ml K r 9 
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where eq.10 corresponds to the temporal coordinate, eq.ll to the radial one 
and eq.12 to the angular coordinates (both produce the same equation). 
Now, from eq.10 we find the coefficients of the A expansion, from eq.ll we 
find the B coefficients and finally, eq.12 is found to be compatible with the 
other two. The values found for the coefficients go as follows:ao = 1, a± = 
2M/M P , a k = af for A; = 2,3,4,5 and a k = a\ + 167rcm^ 4 (54 + (k — 6)5) for 
k = 6, 7, 8, 9 and b — 1, b\ — —a±, b k = for k = 2, 3, 4, 5, 6, b 7 = — 16na5al 
and b & = b 9 = 0. As expected, the first corrections to the Schwarzschild 
metric are of the order of (L p /r) 6 . Note also that now A is not equal to 
B^ 1 and this fact will have important consequences that will be discussed 
later. With regard to the validity of the coefficients obtained above, the 
approximation we have used introducing the Schwarzschild metric in the 
RHS of eq.6 gives correct values for the coefficients up to terms of order 
9. However, when one consider the effect of possible counterterms in the 
lagrangian with eight or more derivatives, they could affect coefficients of 
terms of order seven or higher. 

To study the effect of the new a term on the static isotropic solution one 
could consider an alternative approach by making a perturbative expansion 
in the a coupling constant instead of an expansion in L p /r. To the lowest 
order, it is possible to find: 

1 ,2Af 9/ 1 . R . 54 10M/(rM2) N , 
A (r) = ~, nr + 167ra(— ) 2 (— ) 6 (- + n '\m^ ) (13) 



rMj, ^ ^ rMf, V rMl > 

2M .. 2.1/ 



B(r) = 1 - ^ - 807ra(^) 3 (— — ) 7 (14) 
v 1 rM 2 P K M P J rM P J v ; 

It is easy to see that the behaviour of this solution close to the gravitational 
radius r s = 2M/Mp is now controlled by the a term. For instance, it is 
now possible to find a value near r s where A(ta) = but for the ordinary 
Schwarzschild solution such a point does not exist. In this case the A function 
diverges and the B function vanishes at the gravitational radius of the body 
i.e. at the event horizon A(r s ) = oo and B(r s ) = 0. Similarly, one could 
ask where this happens in our modified solution. Then it is not difficult 
to see that for general a the point r B where B(r B ) = is different from 
the point where A diverges which continues being the gravitational radius 
A{r s ) = oo. For example, for 79 | a | 7i(Mp/M) 4 << 1 (which is a very good 
approximation for macroscopic bodies) we find, for positive a, that rp > 
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r s > r A with r A ~ (2M/M 2 P ))(\ - a57r(M P /M) 4 ) and r B ~ (2M/Mj,))(l + 
o;57r(Mp/M) 4 ) (where A vanishes at r^)- For negative a we have r A > r s > 
r B being and r B given by the same equations. 

More surprisingly, when the scalar curvature is computed for the modified 
metric, it is found to be divergent at r A and r B . This shows the existence of 
an intrinsic singularity at least in the most external radius. In addition, the 
signs of these divergences are different in such a way that in the limit a going 
to zero (where the standard Schwarzschild solution must be found) we see 
that ta and r B go both to r s . In this limit the corresponding two divergences 
of the scalar curvature cancel with each other because of the different sign, 
the singularity disappear as it must be and the A = B^ 1 relation is recovered. 



The modified Robertson- Walker metric: 

The most general isotropic and homogeneous metric can be written as: 

dr 2 = dt 2 - a 2 (t)( , - + r 2 d9 2 + r 2 sin 2 ed<j> 2 ) (15) 
1 — kr z 

where k = —1, 0, 1 corresponding to an open, flat or closed space and a(t) is 
the scale parameter which must be determined by the equations of motion. 
Any isotropic and homogeneous distribution of matter will be described by 
a energy momentum tensor that takes the form of a perfect fluid: 

= (p + p)U"U u + pg^ (16) 

where XJ^ is the velocity four vector, p is the density, and p is the pressure. 
The energy conservation equation reads: 

^ = "3pa 2 (17) 
da 

and finally, the state equation of matter is needed (p = p(p)) to specify 
completely a(t) from the equation of motion and the initial conditions. In 
the case considered in this paper, the equations of motion are: 
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Mp, 3 / n ol , _a" 3 _(a /2 + A;) 3 

' ' a + *)) + 77j(- 24 -T + 12 -« + (18) 



16tt v a 2V " My a 3 a 6 



which is the equation corresponding to the temporal coordinate. The three 
spatial coordinates give rise to only one independent equation. 



Ml,2a" , g» + k , , a rin (a' 2 + kf 

+ ) + W U a° (19) 

, a /w a w + a"' 2 a" (a 12 + k)(ha' 2 + k) a' 2 {a' 2 + k) 2 ^ p_ 
" 72( 3a 2 3a5 + )} + 2 " ° 

Using the energy conservation equation it is not difficult to see that the 
derivative of the first equation above is in fact the addition of the two equa- 
tions. Therefore, we can take eq.18 as the only independent equation. How- 
ever, this equation seems to be very difficult to solve. In order to find some 
exact solutions we have tried to look for a(t) functions which also satisfy the 
relation 

a' 2 + k = Ca n (20) 
The eq.18 forces n to be equal 2 and: 

^C + ^ = P - (21) 
16tt M 2 P 6 y ! 

Thus p must be a constant and this fact, together with the conservation of 
the energy equation imply the state equation p = —p. This means that the 
kind of solutions we are dealing with can describe matter only in the form 
of a cosmological constant p = \Mp/8n. In any case we can use eq.21 to 
obtain the constant C in terms of a and A. Independently of this value, eq.20 
can be easily integrated with n = 2 and the following solutions are obtained: 
a) In the case C = the eq.20 reads a' 2 = —k so for k = a is constant 
and for k = —1, a(t) = ±t + t . This case corresponds to the standard 
Einstein equation with a = in absence of matter. 
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b) For C > eq.20 reads a' 2 = Ca 2 - k. For k — 1, k — -1 and fc = 
the solutions are respectively: 

a(t) = —=cosh(VC(t - t )) (22) 



a(t) = -^=sinh(\ r C(t - t )) 
a(t) = a(t )exp(VC(t - *„)) 
c) For C < there is a solution only in the case k = — 1 where 

a(t) = -±=sin(\/=C(t - t )) (23) 

We would like to stress that all these solutions satisfy exactly eq.18 with the 
Robert son- Walker metric provided that the C is related to a and A through 
eq.21 with p = \Mp/8ir. Therefore now we will write C in terms of the above 
mentioned constants. In particular we have the important case of no matter 
and no cosmological term i.e. pure gravity where we have three solutions: 
C = 0, C = ±^J—Mp/6A7ra. Thus, in this case the C different from zero 
solutions exist only for negative a. 

In the case of general A one must consider the discriminant of the eq.21. 
The sign of this discriminant is in fact determined by a and A. In any case, 
it is always possible to find an analytical formula for the roots of eq.21 in 
terms of a and A. 



Conclussions: 

At low-energies gravitation can be described by an effective lagrangian 
which is an expansion in the number of derivatives over M P . The first term 
of this expansion is the standard Einstein lagrangian but those with higher 
derivatives have undetermined coefficients that finally could be computable 
from the underlying theory of gravitation in a similar way as the coefficients 
of the low-energy effective lagrangian for strong interactions could, at least 
in principle, be obtained from QCD. 

Independently of which are the real values of these coefficients, the con- 
sistence of the theory at low energies requires the aR a ^ 5 W 5 ap R ap ap term 
considered in this work to be present in order to make it renormalizable in 
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the sense of the phenomenological lagrangians. Then it seems to be interest- 
ing to ask about the kind of effects that this term could produce and how it 
can modify the standard predictions of Einstein gravitation. 

In the case of the gravitational field produced by a static and spheri- 
cal body, we obtain that the modifications to the large distance field are 
rather weak and completely unobservable. However, in the region close to 
the Schwarzschild radius, the qualitative behavior of the field changes dra- 
matically since the event horizon is transformed in an extended and spherical 
singularity. In fact, as was already discovered in [11], this is a very generic 
effect and it will appear virtually in any modification of the field equations 
and not only in the one considered in this work. We think that, at least 
to our knowledge, the physical consequences of this fact have not been yet 
properly studied in the literature. 

In the case of the cosmological solution we have shown that inflationary 
solutions exist for the modified equations of motion which do not require the 
somewhat artificial introduction of the inflaton field (this interesting phe- 
nomenon was also found in [12] for the case of an extra R 2 term in the 
gravitational lagrangian). Of course, in the real case one would need to in- 
clude the effect of higher order terms in the equation of motion but our simple 
solutions show that pure gravitation could produce inflation and this could 
be relevant in the very early history of the universe [13]. 

Note added 

When this work was completed, we noticed reference [14] where the au- 
thors study the effect of the six derivative terms on the Hawking temperature 
of black holes. 

Aknowledgments 

This work has been partially supported by the Ministerio de Education 
y Ciencia (Spain) (CICYT AEN90-0034). 



12 



References 

[1] E. Alvarez Rev. of Mod. Phys. 61, 561 (1989) 

[2] Superstring theory , M.B. Green, J.H. Schwarz and E. Witten. 
Cambridge University Press (1987) 

[3] S. Weinberg, Physica 96A 327 (1979) 

[4] J. Gasser and H. Leutwyler, Ann. of Phys. 158 142 (1984) , Nucl. 
Phys. B250 465 and 517 (1985) 

[5] A. Dobado and M.J. Herrero, Phys. Lett. B228 495 (1989) and 
B233 505 (1989) 

J. Donoghue and C. Ramirez, Phys. Lett. B234 361 (1990) 

[6] N. Nakanishi and I. Ojima, Phys. Rev. Lett. 43 91 (1979) 

[7] G. 't Hooft and M. Veltman, Ann. Inst. H. Pomcare2Q 69 (1974) 

[8] M.H. Goroff and A. Sagnotti, Nucl. Phys. B266 709 (1986) 

[9] D.M. Capper, J.J. Dulwich and M. Ramon Medrano, Nucl. Phys. 
B254 737 (1985) 

[10] Gravitation and Cosmology, S. Weinberg, John Wiley & Sons 
(1972) 

[11] W. Israel, Phys. Rev. 164 1776 (1967) 200 

[12] A.A. Starobinsky, Phys. Lett. 91B 99 (1980) 

[13] A. Dobado and A. Lopez, in preparation. 

[14] M. Lu and M.B. Wise, Phys. Rev. D47 3095 (1993) 



13 



